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We use the Saez-BaUester theory on anisotropic Bianchi I cosmological model, with barotropic 
fluid and cosmological constant. We obtain the classical solution by using the Hamilton- Jacobi 
approach. Also the quantum regime is constructed and exact solutions to the Wheeler- De Witt 
equation are found. 
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r—{ ■ Usamos la teoria de Saez-BaUester en el modelo anisotropico Bianchi I con un fluido barotropico 

y constante cosmologica. Obtenemos las soluciones clasicas usando el enfoque de Hamilton-Jacobi. 
El regimen cuantico tambien es construido y soluciones exactas a la ecuacion de Wheeler-DeWitt 
son encontradas. 
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I. INTRODUCTION 
> 

^^ ' Saez and Ballester [l| formulated a scalar-tensor theory of gravitation in which the metric is coupled with a 

CO ' dimensionless scalar field. In this direction, many works in the classical regime have been done [2t5|) yet a study 
C^ ] of the anisotropy behaviour trough the form introduced in the line element has been conected. In this theory the 
f^ , strenght of the coupling between gravity and the scalar field is determined by an arbitrary coupling function w. In 
spite the dimensionless character of the scalar field, an antigravity regime appears, this suggests a possible way to 
solve the missing matter problem in non-fiat FRW cosmologies. However, Armendariz-Picon et al called this scenario 
5_j ■ as K-essence [6[, which is characterized by a scalar field with a non-canonical kinetic energy. Usually K-essence models 
■ - - ' are restricted to the lagrangian density of the form 

\2 



d*xy3^f(^)(v0r, (1) 

one of the motivations to consider this type of lagrangian originates from string theory p]. For more details for 
K-essence applied to dark energy, you can see [g] and reference therein. 

On another front, the quantization program of this theory has not ben constructed, because should be dark to 
build the ADM formalism. Thus, we transform this theory to conventional one, where the dimensionless scalar field 
is obtained from energy-momentum tensor as a exotic matter, and in this sense, we can use this structure for the 
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quantization program, where the ADM formahsm is well known for different classes of matter 

In this work, we use this formulation to obtain classical and quantum exact solutions to anisotropic Bianchi type 
I cosmological model, including a cosmological constant A. The first step is to write Saez-Ballester formalism in 
the usual manner, that is, we calculate the corresponding energy-momentum tensor to the scalar field and give the 
equivalent lagrangian density. Next, we proceed to obtain the corresponding canonical lagrangian Ccan to Bianchi 
type I through the lagrange transformation, we calculate the classical hamiltonian "H, from which we find the Wheeler- 
DeWitt (WDW) equation of the corresponding cosmological model under study. We employ in this work the Misner 
parametrization due that in natural way appear the anisotropy parameters to the scale factors, and we can to analyze 
its behaviour in easy way. 

The more simple generalization to lagrangian density for the Saez-Ballester theory ^1|] with cosmological term, is 

£geo = (R-2A-F((/))0,^,/)''^), (2) 

where (f)'^ — g'^°'(f)^cn R the scalar curvature, F(0) a dimensionless functional scalar field. In classical field theory with 
scalar field, this formalism corresponds to null potencial in the field (p, but the kinetic term is exotic by the factor 
F{<f). 

From the lagrangian ([2]) we can build the complete action 

I = / V-K'Cgco + £mat)d'^x, (3) 

where Cmat is the matter lagrangian, g is the determinant of metric tensor. The field equations for this theory are 

Ga/3 +ga/3A-F(^) f (/)^Q0,/3 - -ga;3(/)_-y0''' j = SttGTq^, 

2F(0)0'r„ + ^0,^0^^ = 0, (4) 

where G is the gravitational constant and as usual the semicolon means a covariant derivative. 

These set of equations Q can be obtained using the equivalent lagrangian as a matter and energy-momentum 
tensor for this field 6, 



L^ = F(</))g"^ 



■ a'P,^: 



Ta0('/') = F(0) f (?!),„(?!),^ - -ga;3(/i,^0''' j , (5) 

in this way, we write the action ([3|) in the usual form 

I = / y^ (R - 2A + /:,„at + l^,i.) d*x, (6) 

and consequently, the classical equivalence between the two theories. We can infer that this correspondence also is 
satisfied in the quantum regimen, because only is modified the hamiltonian constraint [9|. 

This work is arrangede as follow. In section II we construct the hamiltonian density for the cosmological model. 
In section III the classical solutions using the Hamilton- Jacobi formalism are found. Here, we have used a barotropic 
perfect fluid as a matter content and a cosmological constant, obtaining the solutions for differents epoch of the 
evolution for this cosmological model. In Section IV the cuantization scheme, obtaining the corresponding Wheeler- 
DeWitt equation and its solutions for different values for the 7 parameter. Finally, the section V is devoted to 
discussion. 



II. THE HAMILTONIAN DENSITY 

The line element for the cosmological Bianchi type I has the form 

where N is the lapse function, /3± are the corresponding anisotropic parameter in the scale factors, fl play the role 
as the scale factor like to flat Fricdmann- Robertson- Walker cosmological model (e^ = A). The total volume for all 
diagonal Bianchi cosmological models is given by the expression V = e^^^'^^that will appear in the solutions for all 
parameter in this theory. 

Then, the corresponding lagrangian density in this theory is 

wich can be rewritten in the canonical form, 

£ca„ = nof2 + n^+/3+ + n^_/3_ + n^0-N?{, (8) 

with H as the hamiltonian density, and the momentas are defined in the usual way Il^i — 44, where q' = (il, /3+, /3_, (f)) 
are the field coordinates for this system, 

dn ^ 2 N 12 ' ^ ' 

n,, .^.^12^!, ^ ]W_,^^ _, /3,.-Ne-::^, (lo) 

n, ^^^-12^ ^ n^_g^^ ^ /^^-Ne-::^, (u) 
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The matter is introduced as a barotropic perfect fluid P = "fp with 7 a constant between — 1 < 7 < 1, who energy- 
momentum tensor is given by the following expression T^j^ = (p + P)U^U,y — g^i/P where U^ is the four-velocity, p is 
the energy density and P is the thermodynamic pressure in the fluid, respectively. Using the covariant cstructure in 
this tensor, we obtain the differential equation 3flp + 3flP + p = and the solution 

p = p^e-'^^^+'^K 

with pj a constant for the corresponding scenario. 

So, we obtain the following cxprcsion for the canonical lagrangian 



^r„-30 



24 



no - n^^ - nl_ + r^ - iM-^Gp^e'^^^^-^'^ + 48Ae6" 
and 



6n2 

'^+ " ''^- + F(0) 



„-3n / f;rr2 \ 

^ = V T" ""'+""'- + F(^ " 384.0^,6^^(1-'') +48Ae<'^ , (13) 



and using the lagrange equation for the fleld N, we get the Bamiltonian constraint 

6B2 



H = Ill- B2^ - B2_ + —A _ 3847rGM7e3"(i-'') + 48Ae6^ = 0. (14) 



III. CLASSICAL REGIMEN: HAMILTON-JACOBI APPROACH 

Employing the Hamilton- Jacobi formulation, where the momentas are Ilq = -g-^, where Sq is the superpotential 
function, the hamiltonian takes the following form 






dS 



/3+ 



d/3. 



dS 



P- 



a/?- 



6 /dS 



F{(t>) V 90 



3847rGAi^e3"(i-'') + 48Ae^" = 0, 



solving for the variable fl we have 



^ ) - 3847rGM^e3"(i-'') + 48Ae6" = 



fdSp^Y , fdSp_Y 6 fdSA\ 2 



using ([HI, we obtain the following integral equation depending on time parameter Ndt — dr, we have 



dSo 



At = 



dn 



—ii = ./^2 + 3847rGAi^e3f^(i-'r) - 48Ae6^ = 12e'^"-- 



dr 



12 



V^ 



2„-6f^ 



- 3847rGAi^e-3^(i+7) - 48A 



zdr^ 



(15) 



(16) 

(17) 



where ^ is a separation constant. This equation do not have a general solution, but we can solve for the particular 
values of the 7 parameter, which we present in the Table I. 



Case 


n{r) 


Inflation, 7 = — 1 

b-i = 3847rG/i_i - 48A > 


iLn 




Dust, 7 = 

bo = 3847rG/io, A < 


iLn 


1 ^V-SA&T bf, \- ,^2 




(° 4v^mj ^^ 


16\/-3AcV-aA/ix 


Stiff matter, 7 = 1 

bi =^2_^3847rG;Ui, A<0 


iLn 


1 /„v'-3AAt ih,r~^ 


3AAt\ 


_16V-3A {^ ^^^^ 



Table I. Solutions for Q in the scenarios 7 = — 1, 0, 1 



The corresponding solutions for the anisotropic functions and the field (j) appear in the Table II. For the field (j) we 
consider 

6 /9Srf^' 



F(<^) V 9(1^ 



= -^^ + kI + k1^9^ 



6 dr 



with 6 = — ^ + K^ + K_, with the cuadrature solutions 



vm' 



2VG 



/e-3"(r)dr. 



(18) 



Now considering the original Saez-Ballester theory, F(0) — wcf)"'^, the corresponding solution for all m are 



Case 



7 = -1, w > 



fc+ = — v3fc- 



/3±(t) 



T^tanh-i 



2« 






Exp 



2? 



Exp 



^^«=tanh-i 



Exp(^Arj 



, m^-2 



m=-2 



7 = 
w > 



±^tanh-i (^^ 



+4\/^3Xc 



8v^^3A{ 



± 



26l(in+2) 



Cvew 



Exp 



±. 



48 



"fvew 



tanh 
tanh" 



-bo+4\/^3Xe'' 



8\/^3AC 
-bn+4x/^^3Ae^ 



^1 [^ -bn 



-) 



8%/^^3AJ 



, m^ -2 
m=-2, 



7 = 1, w > 



^^tanh-i'° 



"3v1^ 



2VT^ 



T 



29(m+2) 

v^FbTU 



tanh" 



Exp 



46 



%/6br 



=tanh 



^1 /e^^ff^A 
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m7^ -2 
m=-2 



TafeZe //. Solutions for the anisotropic variables (3± and field (/) in the scenarios 7 = —1,0, 1 

These set of solutions satisfy the Einstein field equation (|4]), which were checked using REDUCE package. It is 
interesting to note the behaviour to the cosmological constant in this theory. In the inflationary scenario, we have 
a positive value in such a way that the universe has a bigger growth (but next change to negative one, where the 
universe have a small growth). 



IV. QUANTUM REGIMEN: WHEELER-DEWITT EQUATION 



For quantum regime, we calculate the Wheeler-DeWitt equation 



HM- = <^ n^, - nj^ - III 



SSAnGn^e^"^'^-'^'' + 48Ae^" ^ * = 



where the momenta operators Ilq = —ifi-0-, ^ is the wave function of the universe, also we choose h — 1, thus 



i?* = 



d \ d^ 52 Q^-,n Q2 






- 3847rGAiT,e3"(i-'') + A^ke^^X * = 0, 



where we have used I — ^^ + Q-^ ) for solving the factor ordering problem. Applying the separation method, using 
for the wave function 



* = ^(fi)6(/3+)C(/3_)P(0) 



we obtain 



1 60-" d^V 



A'dn^^ A'dn^B'dM^Cdf]^ V uj ^ 



- 3847rG'/i^e3^(i-^) + 48Ae^" } = 





1 d^B 




BdPl 




1 d^C 




Cdpl 


1 6(1)- 


-"' d^V 



yielding the following set of differential equations 

1 A2r3 

ai (20) 

ai (21) 

•n ^w>2 a^, (22) 

where a| = — a^ + aj + ag, with a| separation constants. The choose of sign for these constants is arbitrary, in absence 
to initial conditions of our universe, studied under this cosmological model. 
The solution for the equations (^01 I^Tj) have the generic form 

g^g±a2/3+^ C = e±^^''", (23) 

The solution for the equation (I22p is more complicated, because depend to the constant m, which is a parameter 
to the Saez-Ballester theory. This equation is rewritten as ^^ — " g T> — 0, which is analog to the equation find in 
the reference 1(/|, y" — ax^y = 0, with a = ^^^ and n — m. 

1. case m=-2 correspond to the Euler equation, who solution have the following structure [lO| 



2? = V^ <( ci + C2Ln0 si a = -| (24) 

Ci sin (/iLn0) + C2 cos {^hnq 



si 


a > 


1 

! 
! 

4 


si 


a — 


si 


a < 



where /i = | \/\l + 4a| > 
2. In the case m=-4, we introduce the following transformation z — ^ and ^ = D, yielding to differential equation 

dz' 



more easy to solve 4-^ — au = 0; so the solution become as 



{Ci sinh (Va^) + C2 cosh (y/Ecf)) si a > 

Ci sin ( Y^|a| 1 + C2 cos ( y^ j si a < ^ ' 

the case a=0, is descarted, because this imply that cu — 0, yielding to the Einstein theory. 

3. When the m parameter satisfy the relation — r^ = 2n + 1, where n is an integer, the general solution take the 
form 



V = 



n+l . . . . 

K I si n > 

6 q ' — 



i^'-'''^) [D6Exp(yff)+D7Exp( 
{^'^''^i^Y'^ [D6Exp(yff ) +D7Exp(-yff ) 



(26) 
si n < 



where Dq, D^ are integration constants and q = ^-^^^ = — — 



2 2n+l ■ 



4. General solution for any m, the solution is expressed in terms of the Bessel function and modified Bessel function, 
for the field 6 



2?= V0Z. (^0^), (27) 



where Z^ is a generic Bessel function, :/ = ^ is the order to the corresponding Bessel function, q = i (m + 2). If 
a < imply that w < 0, Z,y become the modified Bessel function, {I,y,K^). When a > 0, — > w > 0, Zi, — > {Ji,,Yi,). 



On the other hand, the equation ([191) does have not general solution, then this is solved for particular cases in the 
7 parameter, and for this, is rewritten in the following form 



^ - Q^ + (3847rGA.^e3"(i-^) - ASAe'"" ~ a^l) A ^ 0, 
dS2^ diz V / 



1. Any factor ordering Q and the inflation phenomenon! 7 = — 1 

'df) 



d02 



making the transformations z 



A^A A A 

Q— - + [b_ie^" -(il]A = 0, b_i = 3847rG/^_i - 48A 



(28) 



(29) 



— v5lc30 



and A = ze^{z) we arrive at Bessel differential equation for the 



function $. With this the general solution become [ld| 



A = 



±^3n 



Z. 



Igsn 



i. = ±-\/Q2 + 4a2 



(30) 



where Z^ is a generic Bessel function. If b_i > 0, we have the ordinary Bessel function, in other case, will be 
the modified Bessel function. 



2. factor ordering Q=0 and 7 = 



making the transformation z 



dM 



a3Sl 



(48Ae6" - boe^" + aj) A ^ 



and R — z 3 A, is carried to equation 9z^^ + 9 (|ai + l ) 
(48Az — bo) R = who solution is constructed by the degenerate hypergeometric function Fi(a, b; z) 10| 



A 



^aiO 



Exp 



6Ae 



aiO 



312 



, B k 2 e 

Fi — ^,-al + l; — 
' 18k ' 3 A 



(31) 



dR 

dz 



(32) 



where A = -^, k = f^eA, B(k) = 9k(|ai + 1) + bo. 
3. Any factor ordering and stiff matter 7=1 



bi = 3847rGyLti - aj 



(33) 



in similar way that the first case, the transformations z — '^X ~^ ^ ^'^^ -^ — z~f:^{z), the differential Bessel 
function appear for the function $, then we have the general solution 

(34) 



A 



4a --e 



,30 



ZA^\l-'-^e''' 



6 



with A < 0, having the ordinary Bessel function; In the case when the factor ordering become Q=0, we have the 
same Bessel function, but the cosmological constant could be positive or negative, yielding to modified Bessel 
function and ordinary Bessel function, respectively, and the order become imaginary in both cases v = ±z-^^. 



V. CONCLUSIONS 

One equivalent density lagrangian was build in order to apply the quantum regime in the Saez-Ballester theory, 
in where the constant w can be used to have a lorenzian (-1,1,1,1) or seudo-lorenzian (-1,-1,1,1) signature when we 
build the Wheeler-DeWitt equation. The values for this parameter in the classical one is dictated when we apply 
the condition that we must have real functions, which is encoded in the parameter a, equations (|24l25p . In this 
sense, the classical and quantum exact solutions were found for the cosmological Bianchi type I model in the frame 
of Saez-Ballester theory for the scenarius in the 7 parameter {—1, 0, 1}. The presense of the exotic field (f) does not 
delay the anisotropic behaviour in this model. Morcless, the classical behaviour of this field for large value in the 
parameter w, it is similar to the anisotropic parameters /3±. 
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